Tunneling from a correlated 2D electron system transverse to a magnetic field 
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We show that, in a magnetic field parahel to the 2D electron layer, strong electron correlations change the 
rate of tunneling from the layer exponentially. It results in a specific density dependence of the escape 
rate. The mechanism is a dynamical Mossbauer-type recoil, in which the Hall momentum of the tunneling 
electron is partly transferred to the whole electron system, depending on the interrelation between the 
rate of interelectron momentum exchange and the tunneling duration. We also show that, in a certain 
temperature range, magnetic field can enhance rather than suppress the tunneling rate. The effect is due to 
the magnetic field induced energy exchange between the in-plane and out-of-plane motion. Magnetic field 
can also induce switching between intra-well states from which the system tunnels, and a transition from 
tunneling to thermal activation. Explicit results are obtained for a Wigner crystal. They are in qualitative 
and quantitative agreement with the relevant experimental data, with no adjustable parameters. 

PACS numbers: 73.40. Gk, 73.50.-h, 73.20.Dx, 73.50.Jt 



I. INTRODUCTION 

Many properties of low density two-dimensional elec- 
tron systems (2DES) are strongly influenced by electron 
correlations Tunneling is one of the most direct 

tools for revealing these correlations, as has been demon- 
strated for systems which display the quantum Hall effect 
Q . In this paper we show that tunneling can directly re- 
veal correlations in a totally different class of systems, 
the low-density 2DES which are far away from the quan- 
tum Hall regime. This happens for tunneling not into 
the 2DES, but from the 2DES into the volume, and for a 
magnetic field B applied parallel rather than perpendicu- 
lar to the electron layer. The effect may not be described 
in terms of a phenomenological tunneling Haniiltonian: 
it is the tunneling matrix element itself that is sensi- 
tive to the electron correlations. As we show, it depends 
strongly, and very specifically, on electron density, and 
also on temperature and the magnetic field. An expo- 
nentially strong deviation of the tunneling exponent in a 
magnetic field from the predictions of the single-electron 
theory have been observed for a 2DES on helium 
However, these observations remained unexplained. 

The field B parallel to a 2DES couples the out-of-plane 
tunneling motion of an electron to the in-plane motion. 
For an isolated electron, which is separated from the 
continuum by a ID potential barrier U{z), see Fig. 1, 
and is free to move in the plane {x,y), this results in 
an exponential suppression of the rate of tunneling de- 
cay. Indeed, when the electron moves by a distance z 
away from the layer, it acquires the in-plane Hall ve- 
locity vh ~ (e/c)B X z. The corresponding kinetic en- 



a "magnetic barrier" muj'^z' 
decrease of the decay rate. 



/2. This leads to a sharp 




ergy mw|^/2 



:z /2 is subtracted from the energy 



of the out-of-plane tunneling motion (wc = \eB\/mc is 
the cyclotron frequency), or equivalently, there emerges 



FIG. 1. The geometry of tunneling from a correlated 2DES 
transverse to a magnetic field; Up is the characteristic in-plane 
vibration frequency. 

The electron-electron interaction can totally change 
the above picture. If the electron system is spatially cor- 
related, forming a Wigner crystal (WC) or a correlated 
electron liquid, see Fig. |l|, the tunneling electron transfers 
a part of its in-plane Hall momentum to other electrons. 
This decreases the loss of the energy for out-of-plane tun- 
neling motion The mechanism is similar to that of 
the Mossbauer effect where the momentum of a gamma 
quantum is given to the crystal as a whole. However, in 
the present case the dynamics of the interelectron mo- 
mentum exchange is very substantial. The characteristic 
momentum exchange rate is given by the zone-boundary 
plasma frequency ujp, which is related to the electron den- 
sity nhy ujp = (27re^n'^/^/m)-'^/^. If Up exceeds the recip- 
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rocal duration of underbarrier motion in imaginary time 
tJ^ , the WC momentum adiabatically follows that of the 
tunneling electron. As a result, the Hall velocity is the 
same for all electrons, and vh oc 1/iV {N is the 
number of electrons). The effect of the magnetic field on 
tunneling is then fully compensated. For WpTf ~ 1 the 
compensation is only partial. One can say that tunneling 
is accompanied by creation of phonons of the WC, and 
the associated energy goes towards the magnetic barrier. 
However, the barrier turns out to be smaller than for a 
free electron, and the tunneling rate is then exponentially 
larger. Still, for T = it is much smaller than for B = 0. 

We show in this paper, that unexpectedly, in a certain 
temperature range the _B-induced suppression of the rate 
of tunneling from a 2DES may be reversed, and then the 
decay rate exponentially increases with B. This happens 
because thermal energy of the in-plane electron motion 
is transferred by the magnetic field into the energy of 
tunneling motion. Although the effect is generic, as we 
show below, it does not arise in systems where the tun- 
neling rate can be found using the instanton (bounce) 
technique 0|, which is traditionally applied to describe 
tunneling for B — There are two reasons which re- 
quire to modify this technique. First, the magnetic field 
breaks time-reversal symmetry, and therefore, except for 
the case where the Hamiltonian of the system has a spe- 
cial form 1^, there are no escape trajectories in real space 
and imaginary time, and the system comes out from the 
barrier with a finite velocity [|lO|. Second, for 2D systems 
the confining potential U{z) is usually nonparabolic near 
the minimum, and even nonanalytic, with a step in the 
case of heterostructures and, in the case of electrons on 
helium, the singularity of the image potential. Therefore 
such systems are good candidates for observing magnetic 
field enhancement of the tunneling rate. 
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FIG. 2. Energy transfer from thermal in-plane motion into 
out-of-plane tunneling motion induced by the magnetic field 
(schematically) . 

The crossover from suppression to enhancement of tun- 
neling by the field occurs for a temperature Tc which is of 
the order of the reciprocal imaginary duration of under- 
barrier motion Tq~^ for i? = in the ground electron state 
[in what follows we use units where ?i = fcs = 1]. If the 
tunnehng rate for i? = is written as Wq oc exp[— 26*0], 
then To = dSo/dEg, where Eg is the ground-state en- 
ergy of out-of-plane motion in the potential well of U (z). 



The tunneling probability increases if the in plane en- 
ergy Emplane is transferred by the magnetic field into the 
tunneling energy Eg, at least in part, as illustrated in 
Fig. |2| The probability to have an energy iJpianc is 
oc exp(— i5piano/r). Therefore the overall probability, 
which is determined by the product of the two proba- 
bilities, depends on the interrelation between T and tq. 

The time tq also often determines the temperature Ta 
for which there occurs a crossover from tunneling decay 
to decay via activated overbarrier transitions for B — 
pT| , p^ . Therefore Ta and Tc are of the same order of 
magnitude. The interrelation between these tempera- 
tures is determined by the parameters of the system, and 
various interesting situations may occur depending on 
these parameters, as we discuss below. For example, the 
logarithm of the escape rate may increase with B even 
for T > Ta, because in a certain _B-range, the rate of 
tunneling from the ground state exceeds the activation 
rate, even though it is smaller than the activation rate 
for B = 0. Similarly, with increasing B there may oc- 
cur switching from tunneling from the excited intrawell 
states (see Fig. ||) to tunneling from the ground state. 

For T < Tc, on the other hand, the tunneling rate de- 
creases with the increasing B. However, for large enough 
B this increase stops, and escape occurs via thermal ac- 
tivation. 

Explicit results on the effect of electron correlations on 
tunneling will be obtained assuming that electrons form 
a Wigner crystal. Because of strong correlations, over- 
lapping of the wave functions of individual electrons is 
small, and electrons can be "identified" . The problem 
is then reduced to the tunneling of an electron coupled 
to in-plane vibrations of the Wigner crystal. As we will 
see, the results provide a good approximation also for a 
correlated electron liquid. 

In Sec. II we formulate the model. In Sec. Ill we ob- 
tain the general expression for the tunneling rate in the 
WKB approximation, with account taken of the discrete- 
ness of the energy spectrum of electron motion transverse 
to the layer. The result can be understood in terms of the 
tunneling trajectory where the duration of motion trans- 
verse to the layer (in imaginary time) is not fixed, it has 
to be found and depends on temperature and the mag- 
netic field. In Sec. IV we analyze the tunneling exponent, 
including the case of T = 0. In Sec. V we discuss tem- 
perature effects and show the possibility of S- induced en- 
hancement of tunneling and of switching between differ- 
ent regimes of escape from the potential well. In Sec. VI 
explicit results are obtained using the Einstein model of 
the Wigner crystal in which all phonons are assumed to 
have the same frequency. Explicit expressions are ob- 
tained for a triangular and square tunneling barriers. In 
Sec. VII we apply the results to electron tunneling from 
helium surface and provide a detailed comparison with 
the experimental data Sec. VIII contains concluding 
remarks. 
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II. THE MODEL 

A 2D electron system displays strong correlations if the 
ratio r of characteristic Coulomb energy of the electron- 
electron interaction e^(7rn)^/^ to the characteristic ki- 
netic energy is large (here, n is the electron density). 
In degenerate systems the kinetic energy is the Fermi en- 
ergy -Kn/m^ whereas in nondegenerate systems it is the 
thermal energy T . An example of a strongly correlated 
nondegenerate 2DES is electrons on helium. The exper- 
imental data for this system refer to the range F > 20 
1^. A classical transition to a Wigner crystal (WC) was 
observed for F « 130 [^3|jlJ]. Recently much attention 
have attracted also strongly correlated low-density elec- 
tron and hole systems in semiconductors, where there 
have been reached the values of F ~ 40 which are ex- 
pected to be sufficient for Wigner crystallization in a de- 
generate system [0. 

The effect on tunneling of the magnetic field B parallel 
to the electron layer is most pronounced if the tunnel- 
ing length L is long, because the in-planc Hall momen- 
tum due to tunneling mujcL is simply proportional to L. 
Respectively, of utmost interest to us are systems with 
broad and comparatively low barriers. Yet in experimen- 
tal systems the barrier widths are most likely to be less 
than lO'^A. Therefore, in order to somewhat simplify the 
analysis we will assume (although this is not substantial) 
that L is less than the average inter-electron distance 
~ ^-1/2 ^j^jg case, since electrons in a strongly cor- 
related system stay away from each other, the in-plane 
electron dynamics only weakly affects the tunneling po- 
tential jl5|. We will neglect the corresponding coupling 
for B — Q. The major effect on tunneling comes from 
a few nearest neighbors, and the presence or absence of 
long-range order in the 2DES does not affect the tun- 
neling rate. Therefore we will analyze tunneling assum- 
ing that the electron system is a Wigner crystal. As 
we will see, the results will indeed depend on the short- 
wavelength modes of the WC, as expected from the above 
arguments, and therefore we believe that the model pro- 
vides a good approximation even where electrons form a 
correlated fluid. 

In a strongly correlated system, exchange effects are 
not significant, and one can identify the tunneling elec- 
tron. Its out-of-plane motion for i? = is described by 
the Hamiltonian 

H, = ^ + U{z). (1) 

The potential U{z) has a well which is separated by a 
tunneling barrier from the extended states with a quasi- 
continuous spectrum, cf. Fig. |. The well is nonparabolic 
near the minimum, in the general case. The metastable 
intrawell states are quantized. We will consider tempera- 
tures for which nearly all electrons are in the lowest level, 
with energy Eg. 



The magnetic field B parallel to the layer mixes the 
out-of-plane motion of the tunneling electron with the 
in-plane vibrations of the Wigner crystal. The full Hamil- 
tonian is of the form 

H = Ho + Hb+H,, (2) 

with 

Hv^^Yl ["^"^PkjP-ki + mw^jUkjU-kj] (3) 

kj 

and 

Hb = \mLolz^ ~ L^ez7V'i/2 ^[B X pk,],. (4) 

kj 

Here, pkj, Ukj, and LO\^j are the 2D momentum, dis- 
placement, and frequency of the WC phonon of branch 
j {j = 1,2) with a 2D wave vector k. We chose the 
equilibrium in-plane position of the tunneling electron 
to be at the origin. Then its in-plane 2D momentum is 
p = Ar-i/2 J2 pkj for B = 0. 

The interaction Hamiltonian Hb (|^) does not conserve 
the phonon quasi-momentum k. The Hall momentum 
of the tunneling electron is transferred to the WC as a 
whole. The term Hb couples the out-of-plane motion to 
lattice vibrations. The problem of many-electron tunnel- 
ing is thus mapped onto a familiar problem of a particle 
coupled to a bath of harmonic oscillators [^6|j|] , with the 
coupling strength controlled by the magnetic field. The 
distinctions from the standard situation stem from the 
non-parabolicity of the potential well near the minimum 
and from the fact that coupled by Hb are the electron co- 
ordinate z and the in-plane momenia of the lattice. These 
quantities have different symmetry with respect to time 
inversion. In the general case [for example, where the po- 
tential energy of the system has odd-order terms in the 
displacements Ukj], the broken time-reversal symmetry 
requires a special approach to the analysis of tunneling 
pOf . The results discussed below can be appropriately 
generalized using this approach. 

For the model (||), the analysis is simplified by the 
structure of the Hamiltonian (cf. ||^). For vibrations with 
the Hamiltonian H^ (^), one can always make a canonical 
transformation from the canonical coordinates and mo- 
menta Ukj and pkj to the new canonical coordinates and 
momenta pkj and — Ukj, respectively. This transforma- 
tion interchanges the time-reversal symmetry of the in- 
plane dynamical variables, it makes pkj and Ukj even and 
odd in time, respectively. Because Hb is independent of 
Ukj and is linear in pkj , in the new variables it takes on a 
more familiar form of a "potential" coupling, with energy 
which depends on dynamical coordinates only and with 
"restored" time-reversal symmetry. 
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III. THE WKB APPROXIMATION 
A. General formulation 

We will evaluate the escape rate W in the WKB ap- 
proximation. The major emphasis will be placed on the 
tunneling exponent. We will assume that the escape rate 
is much less than the intrawell relaxation rate for relevant 
states, and there is formed thermal distribution over the 
intrawell states of the system. This is not necessarily 
true for 2D systems. Our results can be generalized to 
the case of slow intrawell relaxation, see Sec. VII. 

Electron-phonon interaction under the barrier is 
strong. One should therefore think of escape of the 
coupled electron-phonon system. It results from decay 
of the metastable intrawell states a, with decay rates 
Wa- These rates sharply increase with state energies Ea, 
whereas the Boltzmann intrawell distribution exponen- 
tially decreases with Ea- As a result, there is a com- 
paratively small group of states which mostly contribute 
to the escape (for fast intrawell relaxation, the relataive 
population of these states remains unchanged). This al- 
lows one to characterize escape by a single rate W, 



W^Z-^y iy„exp(-/3^„), 

= Caexp[-2S'„(|f,|in)] |V^a(li. 



(5) 



Here, we introduced a vector ^ — {z, {pkj}) with compo- 
nents which enumerate the z-coordinatc of the tunneling 
electron and the "coordinates" pkj of the phonons, Z is 
the partition function calculated neglecting escape, and 
Ca are the prcfactors in the partial escape rates, they 
will not be discussed in this paper. 

The exponents in the rates Wa are determined |0 by 
the wave functions tpaH) at the turning points on the 
boundary of the classically accessible ranges (^f depend 
on a, see below). It is convenient to evaluate ipai^f) 
in two steps, each of which gives an exponential factor. 
The first factor, exp[— S'a(^f , ^in)], describes the decay of 
the wave function under the barrier. Formally, it relates 
"00 (^f) to ipai^in)- The point is chosen close to the 
well, yet it lies under the barrier, so that Sa can be cal- 
culated in the WKB approximation. The second factor is 
ipai^in) itself. The resulting rate should be independent 

of lin- 

We start with the function Sa , ^in) • To the lowest or- 
der in h, for systems with time-reversal symmetry (which 
we "restored" by the canonical transformation) it is the 
action for a classical underbarrier motion in imaginary 
time T = it with purely imaginary momenta |^ 



Pz = idS/dz, = —idS/dp]ij. 



(6) 



As a function of the imaginary time t, the action S'(^, ^in) 
is given by the integral of the Euclidean Lagrangian Le, 



Sa{i,iin)^ / Lsdr ~ Ea 
Jo 



(7) 



The Lagrangian Le is obtained from the Hamiltonian 
(^ using the Legendre transformation L = pz{dz/dt) — 
Ukj (rfpkj /dt) — H], followed by the transition to imag- 
inary time, which gives 



(8) 



Here, 



and Ly is the Lagrangian of the phonons, Ly — X^kj ^kj, 
with 



L 



1 



1 dpkj dp-kj 



dr 



(10) 



The classical equations of motion in imaginary time 
have the standard form 



d 9Le OLt, 



dr d$, d$, 



= 0. 



(11) 



where overdot means differentiation over r. To calcu- 
late the escape rate, one has to find the trajectory which 
goes from ^(0) ~ to the boundary of the classically 
accessible range at a certain time Tf and calculate the 
action Sa along this trajectory. 

If the potential barrier U{z) is smooth, the wave func- 
tion and its derivatives under the barrier have to match 
the WKB wave function in the classically allowed range 
behind the barrier. The matching occurs at a turning 
point of the classical motion (111]) where the derivatives 
of the both wave functions become equal to zero ||l^, i.e. 
for dSa/dz = 95a/9pkj = 0, i.e. 



i(T/) = 0, Pkj(T/) = 0. 



(12) 



Eq. ([1^) is also the condition of the extremum of Sa with 
respect to the points ^ on the boundary of the classically 
accessible range: the escape rate is determined by the 
minimum of Sa on this boundary. A detailed analysis of 
the behavior of multidimensional tunneling trajectories 
in imaginary time for systems with time-reversal symme- 
try is given in Ref. [^ . 

Time-reversal symmetry of the equations ( |ll|) i n coor- 
dinates (z,pkj), together with the condition (|l2|), show 
that, if the equations of motion are extended beyond t/, 
the system will bounce off the turning point and then 
move under the barrier back to the starting point. The 
section of the trajectory for t > Tf is mirror-symmetrical 
to that for t < Tf, 

z{Tf +t)= z{Tf - t), Pki(T/ + t) = Pki(T/ - t), (13) 
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where < t < tj . As a result, the tunncUng expo- 
nent 2Sa can be calculated along the trajectory ( [Tl| ) that 
reaches the turning point at Tf and returns to the well 
at 2r/. 

The time r/ is determined by the boundary conditions 
( p^ ) and by the initial conditions on the trajectory, which 
are given by ipai^in)- If the intrawell dynamics is semi- 
classical, the dominating contribution to the overall rate 
W (|^) conies from the energies Ea for which the duration 
of the tunneling motion = /3/2 [Q. In the general case 
this is no longer true. 



B. The wave function close to the well 

We are interested in the case where the width of the 
quantum well is much less than the typical width L of 
the tunneling barrier. More generally, we assume that 
for low-lying intrawell states n, the characteristic lengths 
l/7„ of localization in the z-direction are 7"^ L. 
Then, even where the effect of the magnetic field accumu- 
lates under the barrier and the tunneling rate is strongly 
changed, the field may still only weakly perturb the in- 
trawell motion. In this case, inside the well and close 
to it, the out-of-plane electron motion is separated from 
the in-plane vibrations. Respectively, the states of the 
electron-phonon system can be enumerated by n and the 
phonon occupation numbers ntj, i.e. a = (n, {rikj}), and 
the energies are 



Ea ~ En 



(14) 



Usually the interlevel distances — En ^ for 

low-lying levels. 

Because of the separation of motions, we can choose a 
plane z = Zjn under the barrier but close to the well, so 
that for ^ w the wave functions ipaiO ^^'^ semiclassi- 
cal and at the same time can be factored. 



exp 



'S'nkj(Pki) 



(15) 



The action S'„^^ determines the dependence of the wave 
function on the phonon coordinates. 

For ^ = ^iii, Eq. ( p^ ) gives the initial values of the 
dynamical variables ^(0) = and ^(0) on the WKB 



trajectory (|1 1|) . In particular if, for z w Zjn, the poten- 
tial U{z) varies over the distance much bigger than l/7n , 
then 



z{0) 



i(0) - ^ = 
m 



2[C/(zi, 



En 



1/2 



(16) 



and 7„ (16) is independent of the exact position of the 
plane z — Zi^. 

It is convenient to write Sn^^ and pkj in Eq. ( p^ ) in 
the energy-phase representation, using the phonon en- 
ergy and the imaginary time Tkj it takes for a phonon 
to move under the barrier from the boundary (2mekj)^^^ 
of the classically allowed region to the given pkj. With 
the Euclidean Lagrangian of the phonons (p^, we have 

for Pkj = [Pkj]in = Pki(O) 



Sn^j (Pkj(O)) = f drLujir) - e)^jTuj, 



(17) 



and 



Pkj(O) = ekj(2TO£kjy/^ coshwkjTkj, (18) 

1 /2 

Pkj(O) = e_kj (2£kjTOWkj) sinhtJkjTkj 
[ekj is the polarization vector of the mode (k, j)]. 

C. A three-segment optimal trajectory 

To evaluate the escape rate W to logarithmic accu- 
racy, one can, following Feynman's procedure, solve the 
equations of motion (O) for the vibration "coordinates" 
Pkj(T) in terms of z{t) and the initial energies e^j and 
phases WkjTkj- Then, from the boundary condition (|T^), 
one can express r^j in terms of other variables, and then 
perform thermal averaging by summing the escape rate 
over nkj with the Boltzmann weighting factor. Here we 
give an alternative derivation, which provides a better 
insight into the structure of the tunneling trajectory. 

We note that, from Eqs.(|), (0), and (|l^), the 

partial escape rate Wa can be written as Wa oc exp(— s^), 
with 



dr L^,(t) + dT L^{t) + Y,^^ dT L^,{t) - 2EnTf - 2 ekj(T/ + Tkj)- (19) 

I 



(the term 7„Zin in (p^), which is small compared to Sq, ^ 
7„L, is incorporated into the prefactor, see Sec. VII A). 

Eq. (^^ can be interpreted in the following way: 
the tunneling electron in its nth state, accompanied by 



phonons, move under the barrier along a classical tra- 
jectory for the imaginary time 2t/. This motion is de- 
scribed by the Lagrangian Le. Before and after that, 
the phonons are moving on their own, disconnected from 
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the electronic z-motion, for times Tkj and with the La- 
grangian Ly, so that the overaU phonon trajectories make 
closed loops which start and end at turning points. 

In the WKB approximation, the sum over the phonon 
occupation numbers n^j of the weighted partial proba- 
bilities Wa (||) can be replaced by the integral over Skj, 
which should then be evaluated by the steepest descent. 
The corrections due to the discreteness of the values of 
Ekj are small provided uj^jTf <ti Sq. From (p^), the ex- 
tremum of exp(— — (3Ea) with respect to £kj is reached 
for 



(20) 



This expression shows that the duration of the free 
phonon motion is the same for all vibrational modes. 
Moreover, the overall duration of the three-segment op- 
timal trajectory of each vibration is 2(Tkj -I- 77) — [3. 
Examples of the trajectories are shown in Fig. 
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FIG. 3. Optimal trajectories of the tunneling electron z(r) 
and the vibrations of the WC for /3 > 2r/ (a) and /3 < 2r/ (b). 
The numerical data refer to the Einstein model of the Wigner 
crystal, with jTh being the vibrational momentum in the Hall 
direction z x B, in units fi'y /2. The arrows show the direc- 
tion of motion along the optimal trajectory when /3 < 2r/. 
The tunneling potential is of the form (|o|), with dimension- 
less cyclotron frequency ujcTo — 2.0, where to = 2mL/j is the 
imaginary transit time for B — 0. The vibrational frequency 
is ujpTo = 1.0. 

For low temperatures, f3 > 2r/, the direction of time 
along the vibrational trajectory does not change, t^j > 0. 
The corresponding branch of the intrawell vibrational 
wave function ( |l5| ) esc exp(— S'„^^. [pkj (0)]) decays with 
the increasing |pkj(0)| in the classically forbidden region 
|Pkj(0)| > (2TO£kj)"^/^. In contrast, for (3 < 2r/, we 
have Tkj < 0. This shows that the extremum over £kj is 
reached if the intrawell vibrational wave function is an- 
alytically continued from the decaying to the increasing 
branch. 

For Tkj = {(3/2) — Tf < 0, the "free-vibrations" term 
Snt,j is negative, it gives rise to the decrease of the tun- 
neling exponent. This is the formal reason why, for suf- 
ficiently high temperatures, an in-plane magnetic field 



can increase the tunneling rate compared to its i3 = 
value by coupling thermally-excited in-plane vibrations 
to tunneling. 

If the intrawell motion transverse to the layer were 
semiclassical, the sum over the energy levels of this mo- 
tion En in Eq. (||) could be replaced by an integral. The 
extremum of the integrand is reached for Tf — (3/2, T]s^j = 
0. This is the familiar result of the instanton theory, 
in which the whole system moves under the barrier from 
the well to the turning point and back over the imaginary 
time (3 . Clearly, in this case one should not expect the 
tunneling rate to be enhanced by a magnetic field. 

In the case of 2D electron systems, the potential well is 
not parabolic, and each term in the sum over n has to 
be considered separately. Except for narrow parameter 
intervals, the contribution of one of them is dominating, 
and the electron tunnels from the corresponding state. 



IV. THE TUNNELING EXPONENT 

Eqs. (j^) and ( p^ ) allow us to write the escape rate as 
a sum of the escape rates for different intrawell states n. 
To logarithmic accuracy 



W oc max exp(— i?„), R„ 



min 7^„[2;(r)], 

2(t) 



(21) 



The functional TZn [z] is a retarded action functional for a 
ID motion normal to the electron layer. It is determined 
by the functional Sa for the nth state from which the 
the dynamical variables of the in-plane vibrations have 
been eliminated in a standard way | [l6t by solving the 
linear equations of motion ( pT| ) for pkj(T), with account 
taken of the expression for the duration of the phonon 
tunneling motion (EO). This gives 



2r, 



dTi 



+7eccW + (/3-2T/)£;„. 



(22) 



(we count z off from Zjn, i.e. Zin = 0). 

The term TZcc gives the retarded action which results 
from the electron-electron interaction. 



2r, 



dri dT2 



xz(ri)z(r2)xkj(Ti -T2) (23) 
The function Xkj(''") is the phonon Green's function, 

Xkj(T) = flkj CXp[cJkjT] + (fikj + I) CXp[-UJ\^jT] 



^ is the thermal occupation num- 
c [z\ can be also expressed in terms 



(rikj = [exp(/3wkj)-l 
ber). We note that TZ, 
of the correlation function x(t) = (pii (t)p|| (0)) of the 
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in-planc momentum p[| of an electron in the correlated 
2DES, 

7^ee[^] = I ' / dTxdT2 z{ti) z{t2)x{ti " Tj). (24) 

^ Jo ^0 

We expect that not only does this expression apply to a 
Wigner crystal, but it also provides a good approxima- 
tion in the case of a correlated electron liquid; it corre- 
sponds to the lowest-order term in the cumulant expan- 
sion of the appropriate propagator. 

The term TZcc is negative. It means that the electron- 
electron interaction in a correlated 2DES always in- 
creases the tunneling rate in the presence of a magnetic 
field. Moreover, when this term exceeds (mti;^/2) / z^dr, 
the tunneling exponent as a whole decreases with the 
increasing B. 

Two physical phenomena are described by the term 
TZcc- One is the dynamical compensation of the Hall 
momentum of the tunneling electron by the WC as the 
electron moves under the barrier in the z-dircction. The 
other is thermal "preparation" of the Hall momentum 
for the tunneling electron, which is then transformed by 
the magnetic field into the momentum of motion in the 
z-direction. We analyze these effects in the following sub- 
sections. 



A. Zero temperature limit 

It would be natural to think that, since tunneling is 
accompanied by creation of phonons for T = 0, then the 
higher the phonon frequency the lower the tunneling rate. 
In fact just the opposite is true. 

The effect of the electron-electron interaction on tun- 
neling, as characterized by TZcc, depends on the interrela- 
tion between the characteristic phonon frequency ojp and 
the tunneling duration r/. When the tunneling electron 
is "pushed" by the Lorentz force, it exchanges the in- 
plane momentum with other electrons. The parameter 
LiJpTf determines what portion of the momentum goes to 
the crystal as a whole during the tunneling (note that the 
tunneling motion goes in imaginary time, and the quan- 
tity Tf characterizes the time uncertainty rather than the 
actual duration of a real process, see Ref. As men- 

tioned in the Introduction, in the adiabatic limit of large 
tOpTf, all electrons have same in-plane velocity, with an 
accuracy to quantum fluctuations. Therefore the Lorentz 
force produces no acceleration, and no phonons are cre- 
ated during the tunneling. The effect of the magnetic 
field on tunneling should then be eliminated. 

These arguments are confirmed by the analysis of 
Eq. (p3|). If the electron system is rigid enough in the 
plane, so that Lu^jTf :s> 1, the major contribution to TZcc 
comes from T1—T2 '-^ lj^^ <^ Tf. Therefore z(t2) « •^('''i), 
so that in TZcc tbc two terms oc cj^ compensate each other. 



The tunneling occurs as if the electron is disconnected 
from the phonons and does not experience a magnetic 
field. However, a simple analysis shows that the electron 
mass is effectively incremented by a _B-dependent fac- 
tor, and the tunneling exponent i? = i?i is appropriately 
renormalized, 



1/2 



(25) 



R = {m* /my'^ Rb=o 



Here, the sum is limited from below by the condition 



> r 



for a Wigner crystal, the dependence of the 



mass renormalization on the cutoff frequency is logarith- 
mic. The tunneling rate approaches its value for B = Q 
with increasing ojp. On the other hand, the slope of the 
logarithm of the tunneling rate as a function of lOc de- 
pends explicitly on Wp, for ojc ^ i^p- This provides a 
means for measuring LOp. 

For ujpTf ^ 1, only a part of the Hall momentum can 
be taken by the electron crystal. The rest goes into the 
in-plane kinetic energy of the tunneling electron, and ulti- 
mately into creations of WC phonons. However, the ma- 
jor contribution to TZ^c still comes from high-frequency 
phonons. It can be shown from (|2^) that this contribu- 
tion monotonically increases with increasing uj^j . This is 
because the more rigid the electron system is, the more 
effectively it compensates the in-plane Hall momentum. 
An important consequence is that, since high-frequency 
vibrations have small wavelengths, the major effect on 
tunneling comes from short-range order in the electron 
system. 

On the whole, for T — 0, the magnetic-field induced 
term in the tunneling exponent is positive, i.e. the tun- 
neling rate decreases with the magnetic field. This can 
be seen from Eqs. (p^), ( p^ ) by replacing z(ti)z(t2) in 
TZoe with (l/2)[z^(ri) -I- z^Jt2] > z{ti)z{t2) and then in- 
tegrating the function x(ri — T2) = exp[— CL'kj"(Ti — T2)] 
over T2 [for the term z^(ti)] or ri [for z^(r2)]. 

Electron correlations exponentially reduce the effect of 
the magnetic field on the tunneling rate in a magnetic 
field. For specific models, the dependence of the tun- 
neling rate on B and the vibration frequencies will be 
illustrated in Sec. VI, and the results will be compared 
with the experiment. 



B. Small phonon frequencies 

The analysis of the tunneling rate somewhat simpli- 
fies in the case of comparatively high temperatures and 
small phonon frequencies, where the vibrations are clas- 
sical and their frequencies are small compared to the re- 
ciprocal tunneling duration, uj\^j (3 , oj-^jT j ^1. In this 
case 



n, 



,[z] = -2mTujlTjz^ 



7 



/ dTz{T) 

Jo 



(26) 
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[since we chose Zin — and z/ > 0, we have T^cc < in 
Eq. (§)]. 

Eqs. (^), (H^) describe also the tunneUng action 
of a single electron, with the Maxwell distribution of 
the in- plane momentum inside the well exp(--p^/2mr). 
The coupling of the z x B component of the momen- 
tum to the out-of-plane motion gives rise to the term 
—2puJcJf^^dTz{T) in the tunneling action [cf. Eqs. (|^), 
(^]. The extreme value of the sum of this term and 
—p^/2mT is just equal to —TZcc[z] as given by (|2^). 

The single-electron form of the tunneling exponent is 
to be expected in the limit of small Wkj , because the dis- 
tribution over in-plane momenta of electrons forming a 
Wigner crystal is Maxwellian, in the classical limit. For 
small uJkjTf the momenta do not change over the tunnel- 
ing duration, therefore only the momentum of the tun- 
neling electron itself is important. The above derivation 
provides an independent test of the derivation used to 
obtain the general expression (|2^), (|2^). 

We note that the action TZoc[z] ( p6^ is still retarded, 
it does not correspond to a local in time Lagrangian. 
The functional form of TZcc remains the same even for 
temperatures T < Wkj provided the phonon frequencies 
are small compared to tJ^ and lOc- In this case T in 
Eq. ( p6| ) has to be replaced by (4A^)^^ X]'^ki(2fikj + !)• 
This factor explicitly depends on the phonon disper- 
sion law, but again, the major contribution comes from 
short-wavelength high-frequency phonons, which are de- 
termined by the short-range order in the electron system. 



V. ENHANCEMENT OF TUNNELING BY A 
MAGNETIC FIELD 



In this section we describe a new effect, the enhance- 
ment of the tunneling rate by a magnetic field parallel 
to the electron layer. Qualitatively, the enhancement is 
due to transferring the energy of thermal in-plane motion 
into the energy of out-of-plane tunneling. On the formal 
level it is a consequence of the increase, with increasing 
temperature, of the absolute value of the term TZcc ( |2^ ) 
in the tunneling action. Since this term gives a negative 
contribution to the tunneling exponent R, the whole in- 
dependent term in TZ becomes negative starting with a 
certain crossover temperature Tc, and then the tunneling 
rate increases with B. The range boundaries where the 
overall escape rate increases with B are not universal and 
depend on the potential U (z) and the phonon spectrum. 
The enhancement occurs in a limited temperature range, 
and may start from B = or have a finite threshold in B. 
However, very strong fields suppress rather than enhance 
escape. 



A. The crossover temperature 

The lower temperature bound of the enhancement do- 
main is the crossover temperature T^. It can be deter- 
mined from the small-_B expansion of the tunneling ex- 
ponent for the ground state [n = g in Eq. , 



RgiuJc) « Rg{0) + Ag(T)L0l, W,To < 1 



(27) 



where tq is the tunneling time in the ground state for 
B = {). The role of the ground state is special in that 
the barrier width is bigger for the ground state energy 
than for the energies of the excited states. Therefore the 
effect of the magnetic field, which accumulates under the 
barrier, is most pronounced in the ground state. 

The value of Ag is given by the terms oc lo^ in the 
action Rg (|2^) calculated along the tunneling trajectory 



zo(t) for B 

lows that Ag > Q for T 
is given by 



0. From the analysis in Sec. IV A it fol- 
0. The crossover temperature 



Ag(T,) = 0. (28) 

For T > Tc the tunneling exponent Rg decreases and the 
tunneling rate increases with B, for small B. 

In the limit of low phonon frequencies, Wkj ^ ^/tq^Tc, 
from Eqs. (|||), (|§) it follows that f3c = l/T^ = 
2toZq^/zq, where zq is the average coordinate z (p6|) for 



the B — trajectory with energy Eg, and Zq is the mean 
square value of z on the same trajectory. 



drzlir) [E^Eg). 



Clearly, in this case (5c < 2to. It follows from Eq. (|22 
that 2ro is also the limiting value of (3c in the opposite 
case of high phonon frequencies, Wkj ^ I/tq. On the 
whole, we have the bounds on temperature for the tun- 
neling enhancement in the ground intrawell state 

(29) 



2to= < (3c< 2to. 



As noted above, zq is nonzero, and generally zq^/zq ^ 1- 
It follows from the above arguments that the value of 
the crossover temperature Tc = 1/Pc decreases with in- 
creasing phonon frequencies, that is the crossover is de- 
termined by high-frequency phonons which, in the case 
of 2D electron systems, have large wave numbers and are 
determined by the short-range order. 



B. Upper temperature limit 

A thresholdless tunneling enhancement starting from 
B — 0, occurs for temperatures bounded from above by 
the condition that the system tunnels from the ground 
state rather than from excited intrawell states or via ther- 
mal activation over the barrier. In principle, even for ex- 
cited states, the tunneling rate may increase with B, but 
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this does not happen for simple model potentials inves- 
tigated below. 

If the tunneling is enhanced only in the ground state, 
the upper temperature bound is often the temperature 
Ti^2 where the probability of tunneling from the first 
excited state, weighted with the occupation factor, ex- 
ceeds that from the ground state, for B = Q. It can be 
estimated for smooth tunneling barriers, where the tun- 
neling duration to{E) for B — often decreases with the 
increasing energy E. In fact, the function tq(E) may 
be nonmonotonic even for simple potentials U{z); a de- 
tailed analysis of this function lies outside the scope of 
this paper, but generalization of the results to appropri- 
ate cases is straightforward. From (p2[), for decreasing 
To{E), switching from tunneling from the ground state 
{n = 1) to that from the first excited state (n — 2) oc- 
curs for the reciprocal temperature 

This value lies between 2ro(-E2) and 2to{Ei). Depend- 
ing on the tunneling potential, /?i^2 can be smaller or 
larger than /3c • If a magnetic field does not increase 
the rate of tunneling from the state n = 2, thresholdless 
tunneling enhancement occurs for Tc < Ti^2- 

Alternatively, for B = the system may switch to acti- 
vated escape over the barrier with increasing temperature 
for T = Ta < Tc- The thresholdless tunneling enhance- 
ment by the magnetic field does not occur in this case. 
However, in this as well as in the previous case there may 
still occur a B-induced enhancement of the escape rate, 
starting with some nonzero B. We note that in the above 
arguments, it was assumed that thermalization inside the 
well is faster than electron escapes. 

C. Field-induced switching between the levels and 
from activation to tunneling 

Even in the temperature range T > Ti^2 a sufficiently 
strong magnetic field can increase the tunneling rate, 
provided T > Tc- This happens if the tunneling expo- 
nent for the ground state Rg{ujc) = Rn=i{'-Oc) exceeds 
that in the first excited state Rn=2{i^c) and its zero- field 
value -Rra=2(0). In a certain temperature range where 
T > Ti^2, the tunneling rate for B = is determined 
by tunneling from the excited state n — 2. This rate de- 
creases with increasing B (the tunneling exponent Rn=2 
increases with B). For some B the exponents i?„=2(^c) 
and Rn=i{cOc) become equal to each other. For larger B 
the system tunnels from the ground state, and the tun- 
neling rate increases with B. 

Similarly, since the activation rate is only weakly af- 
fected hy B, in a certain temperature range where es- 
cape already occurs via activation for B = 0, starting 



with some B it may again go through tunneling from the 
ground state. This happens if the tunneling rate for the 
ground state becomes bigger than the activation rate and 
only happens in a limited range of B, see Sec. VI. For a 
special model the switching is illustrated in Fig. || below. 



VI. TUNNELING ENHANCEMENT FOR THE 
EINSTEIN MODEL OF A WIGNER CRYSTAL 

In what follows we will illustrate the general results 
and apply them to specific 2D systems assuming that all 
vibrational modes have the same frequency, Wkj = iOp, 
i.e. using the Einstein model of the Wigner crystal. 
This is motivated by the fact that the tunneling is deter- 
mined primarily by short-wavelength vibrations, which 
have a comparatively weak dispersion. Correspondingly, 
when we discuss the experiment, we will set ujp equal 
to the characteristic short-wavelength plasma frequency 
[2-Ke^r?!'^ ImYl"^ ^ where n is the electron density. 

A. Triangular barrier 

For electrons above helium surface and in certain types 
of semiconductor heterostructures, the potential Viz) in 
the barrier region (z > 0) is determined by the electric 
field which pulls electrons away from the intrawell states. 
To a good approximation U{z) is then linear in z, 

f/(.) = |-(l-i) (.>0). (30) 

Here, 7 = 71 is the decrement of the ground-state wave 
function "i/'g = "01 near the well, dlnipi/dz — —7 for 
z — 0, cf. the discussion before Eq. ([l6|). The addi- 
tive constant in U (z) is chosen so that the energy of the 
ground state Eg — 0. Then L is the tunneling length 
in the ground state for B = 0. It is determined by the 
pulling electric field. We assume that 7L ^ 1. 

The approximation (^o|) applies only within the bar- 
rier region, where U is determined by the pulling electric 
field, and not inside the well, where U{z) is singular. 
Moreover, it holds provided the width of the tunneling 
barrier is small compared to the in-plane interelectron 
distance n-^/^ [cf. Eq. (||) below]. 

In order to calculate the ground-state tunneling ex- 
ponent, it is convenient to solve directly the equations 
of motion JTT]) with the boundary conditions (p^), (p^), 
( p^ ) , and (|20|). For a triangular potential, these equations 
are linear. This allowed us to obtain for the tunneling ex- 
ponent a simple expression 

R = -l^pT-rod + 3fpTrod(l - Trod) COth[tJp/3/2 - VpTrcd] 

+ 3 + 3Trcd{l^^ - 1), Rg = 2jLR/3l^^. (31) 
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Here, i^p — ujpTQ and i^c = l^cTq are the dimcnsionless in- 
plane and cyclotron frequencies scaled by the tunneling 
duration tq for B — 0, and i/^ — + v^. 

The quantity Tjod = 'Tf/'To in Eq. (^) is the reduced 
tunneling duration. It is given by the equation 

[(1 - Trcd)VpV^ COth[tJp/3/2 - UpTrcd] - ^l] tanh i/Trcd 



r 2 2 



(32) 



In the limit T ^ 0, Eqs. (|l|), (|32h go over into the re- 
sult obtained earlier (in Ref. [ |1C| we used loq and I'o 
instead of and Vp, respectively). In this limit, the role 
of the many-electron effects is particularly important. In 
the single-electron approximation {ujp = 0) the tunneling 
duration Tf and the tunneling exponent Rg diverge for 
^ Tf^^ 1^- This happens because the effective single- 
electron potential U{z) + {l/2)muj'^z'^, which takes into 
account the parabolic magnetic barrier, does not have 
classically allowed extended states with energy Eg = 
behind the barrier. 



o 
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FIG. 4. The tunneling exponent in the ground state for a 
triangular potential barrier (^o|) as a function of the phonon 
frequency uip in the Einstein model of the Wigner crystal for 
WcTo = 2. The time tq — mL/'y is the duration of tunneling 
for B — and T = 0. The curves 1 to 3 refer to reciprocal 
temperatures /3/ro = 7, 5, 3. The dashed line is the result of 
the direct variational method, with one variational parameter 

Tf- 

The inter-electron momentum exchange makes tunnel- 
ing possible for all B. For lUcTq > 1 and T — 0, the tun- 
neling exponent is a steep function of the exchange rate 
ujp in the limit of slow exchange, ujpTQ <C 1. In the oppo- 
site limit of the fast momentum exchange, Up ^ Tq^,lUc, 
from Eqs. (pl|), (^) we obtain Tjod ~ 1 [he., tj « tq], 
and Rg « A-yL/S. These are the values for tunneling for 
B — 0. The overall dependence of the tunneling expo- 
nent on ujp for T = is shown in the inset of Fig. ^. 

For a given magnetic field, the dependence of the tun- 
neling exponent Rg on the frequency LOp becomes much 
less steep with increasing temperature, as seen from 
Fig. ^ For large WpTQ, LUp/S, the curves for different 



temperatures merge together and approach the B = 
asymptote. 

The value of Rg can be calculated independently 
from the functional TZn (|2^) using the direct variational 
method. Even a simple approximation where z(r) is 
quadratic in t, with the only variational parameter being 
the tunneling duration r/ , gives a reasonably good result, 
which is shown in Fig. by a dashed line for (3 = 3tq. 
Such calculation gives a good approximation for higher 
temperatures, and also for lower temperatures but not 
too small ujqTo. For low temperatures and small luqTq the 
trajectory z(r) is strongly nonparabolic, and more then 
one parameter is required in the variational calculation. 



Field-induced tunneling enhancement and 
switching to tunneling from activation 



The explicit expression for the tunneling exponent ( |3l| ) 
allows us to analyze the effects of tunneling enhancement 
and magnetic field induced switching to tunneling, which 
were discussed in Sec. V. In the small-B limit, where 

LUc < LUp,T, 



,iQ^, the tunneling exponent Rg{B) is seen 
from Eq. (^TI) to be quadratic in B. The coefficient Ag 
in Eq. ( p7[ ) can be easily calculated. From the condition 
Ag ~ Q we obtain the value of the reciprocal temperature 
(3c which corresponds to the crossover from decrease to 
increase of the tunneling rate due to a magnetic field, 



Pc = 2ro H tanh 



Vp[iVp 



(3 -f i^p) tanhz^p 



3t'r, -|- 3 tanh 



(33) 



In agreement with (|29|), (3c monotonically increases with 
Ljp from 5to/3 at Wp = to 2ro for Wp — >• oo. 

The dependence of the tunneling exponent ( plj ) on 
the magnetic field for different temperatures is shown in 
Fig. 1^. Above the crossover temperature {(3 < (3c), R{B) 
decreases, and R{0) — R(B) and the tunneling probability 
increase with the increasing field, for small B. The slope 
dR/BdB oc (3 — Pc for B ^ 0. However, for strong fields 
the tunneling rate decreases with the increasing B, be- 
cause the Hall momentum can no longer be compensated 
by thermal fluctuations. For small ujpTo this happens 
when the typical Hall momentum mujcL becomes com- 
parable to the thermal momentum (2mr)^/^ multiplied 
by the small factor {-/L)~-^/^. This factor comes from 
the fact that the optimal value of the transferred in-plane 
momentum in the single-electron approximation is deter- 
mined by the maximum of the sum 2Se[p] + ip'^/'^.mT), 
where Se[p\ oc 7^ is the single-electron action in the 
magnetic field for given in-plane momentum p. There- 
fore the thermal momentum is scaled by {-/Ly^-^. 
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FIG. 5. The dependence of the tunnehng exponent 
R{B) = Rg{B) on the magnetic field (jsi]) for tJpTo = 1/3 
near the crossover temperature /3c ~ 1.67ro (js^). The curves 
1 to 3 correspond to (/3 — /3c)/ro = 0.2, 0, —0.3 

It is clear from the data in Fig. |^ that, for the bar- 
rier chosen, the magnetic field induced increase of the 
tunneling exponent R is numerically small. However, for 
typical i? > 50 it can still be noticeable, although strictly 
speaking it is on the border of applicability of the approx- 
imation in which only the exponent is taken into account. 
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FIG. 6. Magnetic field induced switching from activation 

(a) and from tunneling from the excited state (b) to tunnel- 
ing from the ground state, for LUpTo — 1/3. In (a), there is 
only one intrawell state in the potential well U{z), and the 
transition to activation for 5 = occurs for P/tq = 4/3. 
The curves 1, 2 correspond to (/3 — I3c)/to = —0.35, —0.4. In 

(b) , the position E2 of the excited level (n = 2) is chosen 
at 0.2'y^ /2m below the barrier top [Ei = 0). The temper- 
ature is chosen at {l3 — f3c)/To = —0.16, so that for 5 = 
the system tunnels from the excited state. The observable 
(smaller) tunneling exponents for a given B are shown with 
bold lines, whereas dashed lines show the bigger exponents, 
which correspond to smaller tunneling rates. 

The expression ( ^l| ) gives the tunneling exponent only 
for low enough temperatures where the system escapes 
from the ground state. For higher temperatures, one 
should take into account the possibility of escape from 
excited states and via an activated transition over a po- 
tential barrier. The positions of the excited levels depend 
not only on the barrier shape, but also on the shape of 



the potential U{z) inside the well. The analysis for a re- 
alistic system, electrons on the surface of liquid helium, 
is done in the next section. Here, in order to illustrate 
different options, we discuss two cases: a narrow well, 
in which case the ground state is essentially the only in- 
trawell state, and a well with a comparatively shallow 
excited state. We assume that the intrawell relaxation 
rate is higher than the escape rate. 

We start with the case of one bound state in the poten- 
tial well. For B = and a triangular barrier U{z) (|3^), 
switching from tunneling to activation occurs here for the 
temperature Ta = l/f3a — (4to/3)^^. This temperature 
is higher than the crossover temperature l/(3c (^3|), and 
therefore there is a region where the enhancement of tun- 
neling by a magnetic field can be observed, as discussed 
above (cf. Fig. H). However, even though for T > Ta 
the B = 0-escape occurs via overbarrier transitions, the 
increase of the tunneling rate with the increasing B can 
make tunneling more probable for sufhciently strong B. 
If the activation rate is independent of B, the overall de- 
pendence of the exponent of the escape rate R{B) on B 
is shown in Fig. In this case, R{0) = 7^/2mT is the 
barrier height over temperature. Switching to tunnel- 
ing and to the increase of the escape rate with B occurs 
where the tunneling exponent Rg{B) as given by Eq. ( |3l| ) 
becomes less than i?(0). 

A similar switching occurs in the temperature range 
where tunneling from the first excited level is more prob- 
able than from the ground state, for 5 = 0. Since with 
increasing B the tunneling rate in the ground state in- 
creases, the system switches to tunneling from the ground 
state starting with a certain value of B. This is illustrated 
in Fig. 

In narrow-well potentials, a magnetic field may 
strongly affect the wave functions with energies close to 
the barrier top. As a result, new bound metastable states 
may appear in a strong field. The field also shifts the 
energy levels of the existing states. The rate of inter- 
level transitions may also change, since the field mixes 
together the in-plane and out-of-plane motions. The re- 
lated effects may become important with increasing tem- 
perature. 



Square barrier: field-induced crossover to 
thermal activation 



In many physically interesting systems, the tunneling 
barrier U{z) is nearly rectangular. This is often the case 
for semiconductor heterostructures, where the barrier is 
formed by the insulating layer. If we count U off from 
the intrawell energy level Eg and set the boundaries at 
z — Q and z — L, the barrier has the form 



U{z) = 7V2m, < z < L. 



(34) 
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Here, I/7 is the decay length under the barrier, cf. 
Eq. (p^), and we have neglected the lowering of the bar- 
rier due to the electrostatic field from other electrons 
at their lattice sites, which is a good approximation for 
nL^ < 1. 

We assume that, behind the barrier {z > L), an elec- 
tron can move semiclassically with all energies. Then the 
decaying underbarrier wave function has to be matched 
to an appropriate propagating wave behind the barrier 
at z = L. In contrast to the case of a smooth barrier, 
because the potential U{z) is discontinuous at z = L, 
the z-component of the momentum should not be the 
same on the opposite sides of the boundary. However, 
the in-plane "momentum" components Ukj, which are 
imaginary under the barrier, still have to be continuous. 
Respectively, the boundary conditions for the tun- 
neling trajectory should be changed to 



z{Tf)=L, Ukj(T/)=0. 



(35) 



In fact, the condition Ukj 



of Pkj for which the wave function is maximal for z 



gives the in-plane values 
L. 

With the boundary conditions (|3^), elimination of 
phonon variables from the Euclidean action Se in the 
tunneling exponent is similar to what was done for a 
smooth barrier. The resulting expression for the re- 
tarded functional TZn[z] coincides with Eq. (p^), provided 
z{Tf + x) is defined as z(t/ — x), for < x < Tf. 

An important feature of a rectangular tunneling bar- 
rier is that the tunneling time to{E) = —dSo/dE for 
B — monotonically increases with energy E. There- 
fore the maximum of the function —(3E — 2So{E), which 
gives the probability of tunneling with energy E with ac- 
count taken of the occupation factor, corresponds either 
to the transition from the ground state or to activation 
over the barrier. Switching to activation occurs for the 
temperature Ta = 7^/4mS'o(£'g) = ^/4mL = (4ro)"^ It 
is lower than the temperature Tc of the crossover from B- 
suppressed to B-enhanced tunneling as given by Eq. (29), 
and therefore we do not expect the crossover to occur in 
systems with a square barrier. 

If the temperature T < Ta, escape for S = occurs 
via tunneling, and its probability decreases with the in- 



creasing B. Starting with some B, where the tunneling 
exponent becomes bigger than the activation exponent 
7^/2mT, it becomes more probable to escape by acti- 
vated transition than by tunneling. To a good approxi- 
mation, the escape rate becomes than independent of the 
magnetic field. 
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FIG. 7. The logarithm of the escape rate R{B) for the 
square potential barrier, cUp = (2ro)~^ (tq = mL/7). 
The value -R(O) is given by the tunneling exponent, 
R{0) = 2So{Eg) = 27L. Curves 1 to 3 correspond to 
/3/ro = 4.5, 5.5, 6.5. The sections of the curves where R{B) in- 
creases correspond to tunneling and are described by Eq. (pq), 
whereas the horizontal sections of the curves correspond to 
thermal activation. Inset: the magnetic field dependence of 
the switching temperature TaTQ. 

The S-dependence of the escape rate for different tem- 
peratures is illustrated in Fig. ^ The results refer to the 
Einstein model of the Wigner crystal. In this model the 
tunneling exponent can be obtained directly from the 
[linear, in this case] equations of motion ( pd| ) with the 
boundary conditions ([l^), (p^), (p5|). It has the form 



Rg = 7L [1-f Trod + J^cK(Trcd)] , 



(36) 



where the function K(Tred) and the reduced tunneling 
time Tred — Tf/ro are given by the equations 



K(Trcd) 



:/c(coshz/Ti.od - 1) 



+ coshi/Ti-cd + tyvp coth[u!pl3/2 — VpTicd] sinhz/Trod 

1 ~ 2coshi/Trd -l- z^Tid sinhi/Trd) — i^^{Trd — l)sinhi'Ti.c 

fct'p (1 - cosh i^Trd) (1 - i^c/'^p) + ^'^rd sinh i/r^d 



(37) 



with i^p — ujpTo, Vc — uJcTq, and v — {v^ + v^Yl"^ . 

The temperature of switching to activation is given 
by the equation = -{^ jlraRg. From (|6|), (|3^), both 
Rg and the reduced time Trcd increase with the mag- 



netic field, and therefore the switching temperature Ta 
decreases with B. However, it follows from the analy- 
sis of the above equations that Ta remains lower than 
l/(2r/), and r/ > To for all temperatures where the sys- 
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tern escapes via tunneling, in contrast to the case of the 
triangular barrier discussed earlier. 

The effect of saturation of the escape rate with increas- 
ing B is not limited to square barriers, of course. For 
strong enough B and nonzero temperatures, the tunnel- 
ing rate becomes less than the activation rate, and the 
system switches to activation; the switching may go in 
steps with increasing B, via tunneling from excited in- 
trawell states. 



VII. COMPARISON WITH THE 
EXPERIMENTAL DATA ON TUNNELING FROM 
HELIUM SURFACE 

Tunneling from a strongly correlated 2DES has been 
investigated in much detail for electrons on helium sur- 
face [p0|j5|. In this system, a good agreement has been 
reached between theory and experiment in the absence of 
the magnetic field, where the primary role of the electron 
correlations is to change the effective single-electron tun- 
neling barrier (see below). As mentioned before, there 
were also done interesting experiments on tunneling in 
a magnetic field. However, the observed strong field de- 
pendence of the tunneling rate dramatically differed from 
the predictions of the single-electron theory and has re- 
mained unexplained 

Electrons on helium surface are localized in a ID po- 
tential box. One side of this box is the image potential 
-A/z, where A = e^ie - l)/4(e-f 1) (e « 1.057 is the di- 
electric constant of helium) , and z is the direction normal 
to the surface. The other side is a high barrier ~ 1 eV 
on the surface (z = 0), which prevents electrons from 
penetrating into the helium. The intrawell states can be 
made metastable by applying a field £± which pulls the 
electrons away from the surface. This field is determined 
by the helium cell geometry and depends on the applied 
voltage and the electron density n, cf. j2^. The overall 
electron potential has the form 

U{z) = -Az^i - \e£^\z - mup-z^ {z > 0). (38) 

The term oc describes the Coulomb field created by 
other electrons at their in-plane lattice sites (the "corre- 
lation hole" [|l^,^). Only the lowest-order term in the 
ratio of the tunneling length L to the interelectron dis- 
tance ri~^/^ has been kept in Eq. (|3^), and we used the 
interrelation 



2m^ ' 



-3 



1/2 



2N ^ 

ki 



^kj 



1/2 



(39) 



for the sum over lattice sites R;. For a triangular 
lattice, uj w (4.456^^3/2/771)^/^ g. The conditions 
1/7 ^ L ^ riT^I'^ are typically very well satisfied 
in experiment, with the decay length I/7 = l/Arn w 



0.7 X lO"*' cm, L - \EgleEx_\ « 7V2m|e£±| - 10"^ cm, 
and •nT^I'^ ~ 10^^ cm [in the estimate of L we used that 
\Eg\ > |e£:j_|/7, 7^^2/72, and that |e£:_L|/7 > w]. 



To compare the predicted dynamical effect of the 
electron-electron interaction with the experimental data 
on tunneling in the magnetic field we use the Ein- 
stein model of the WC. In the equations of motion ( pT| ) 
we assume that all phonon frequencies Wkj are the same 
and set them equal to the characteristic plasma frequency 
— {^■Ke^Jt'l'^ jraYl'^ [yet we use Eq. (39) for the pa- 



rameter a) of the potential barrier ?7(z)]. The numerical 
results change only slightly when the phonon frequency 
is varied within reasonable limits, e.g., is replaced by uj. 



The magnetic field dependence of the tunneling rate 
for the parameters used in the experiment is calculated 
from Eqs. (|ll|), ([l|), (|l|), <^ and is shown in Fig. |. 
The data refer to the values of T where escape occurs 
via tunneling from the ground state. The actual calcu- 
lation is largely simplified by the fact that, deep under 
the barrier, the image potential — A/z in (^) can be 
neglected. The equations of motion (|l^) become then 
linear, and the tunneling exponent i? can be obtained in 
an explicit, although cumbersome form, which was used 
in Fig. ||. The correction to K from the image potential 
is ~ Xj^L^ which is the small parameter of the theory. 
Moreover, since this correction comes from the range of 
small z, where the effect of the magnetic field is small, it 
is largely compensated where i?(i?) — -R(O) is calculated. 
This and other corrections ~ 1/7^ result in changes of 
the theoretical curves that are smaller than the uncer- 
tainty in R{B) — R{0) due to the uncertainties in n and 
£± in the experiment M. 



As seen from Fig. || the dynamical many-electron the- 
ory is in good qualitative and quantitative agreement 
with the experiment, without any adjustable parameters. 
At low temperatures (T = 0.04 K), the many-electron 
tunneling rate is bigger than the single-electron estimate 
[§ by a factor of 10^ for B = 0.25 T. For this tempera- 
ture, the tunneling rate is well described by the T 
limit [ p9[ . The i3-dependence of the tunneling rate is very 
sensitive to temperature. It becomes less pronounced for 
higher T, and the role of dynamical many-electron effects 
becomes less important, too. Interestingly, the theoreti- 
cal data on the ratio of W{B)/W{0) become less sensi- 
tive to the experimental uncertainties in the cell geome- 
try (which determines £±) and the electron density n for 
intermediate temperatures T ^ 0.14 K. This is because 
the corresponding errors in W{B) and W{0) compensate 
each other for such temperatures. 
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FIG. 8. The relative rate of electron tunneling from the he- 
lium surface W (B) /W (0) as a function of the magnetic field 
B for the electron density n = 0.8 x 10* cm~^ and the cal- 
culated pulling field £± — 24.7 V/cm (solid curve). Solid 
lines show how the theory compares to the experimental data 
points of Ref. . The error bars show the uncertainty in the 
theoretical values due to the uncertainty in the parameters of 
the experiment. 

The crossover to magnetic-field enhanced tunneling oc- 
curs for temperature Tc ~ 0.19-?^, for the parameters 
in Fig. ^. The expected increase of the tunneling rate 
with _B for T > Tc is shown in Fig. ||. It has indeed 
been observed in the experiment ||] . The analysis of the 
experiment requires to establish whether, for tempera- 
tures of interest, escape actually occurs via tunneling. 
To that end we note first that, as it follows from a di- 
rect variational calculation, the potential U{z) (|3^), with 
the parameter values specified in Fig. ^, has only one 
metastable intrawell state. 

If the intrawell relaxation were fast enough, the tem- 
perature of the crossover from tunneling to activation Ta 
for B = would be given by the condition that the tun- 
neling exponent Rg be equal to the activation exponent 
(^^max ~ Eg)/T [here, C/max is the maximal value of the 
potential U{z)]. This would give Ta « 0.15 K. However, 
activated escape requires that the in-plane thermal en- 
ergy of an electron be transformed into the energy of its 
out-of-plane motion. This involves a large transfer of the 
in-plane momentum ^ [2m{Uma.->!_—Eg)]^^^ ■ The electron- 
electron interaction does not give rise to such a transfer 
in a strongly correlated system, since the reciprocal in- 
terelectron distance is n^/^ <C [2m(C/max — Eg)Y/^. 

The major process which gives rise to the momentum 
transfer is scattering by capillary waves on the helium 
surface, ripplons Q. Electron-ripplon coupling is weak. 
As a result, the prefactor in the activation rate, which is 
quadratic in the coupling constant, is small. For B ^ Q 
it is ~ j'^T'^/ha 1^^, where a is the surface tension of 
liquid helium. For temperatures T < 0.25 K this pref- 
actor is less than the prefactor in the tunneling rate 
(7i7^/m) exp(— 2) by a factor < 10~^. Therefore the 



crossover from tunneling to activation occurs for higher 
temperatures than it would follow from the condition of 
equal tunneling and activation exponents. 

For the parameters in Fig. ||, the rates of activation and 
tunneling escape become equal for temperatures slightly 
higher than 0.26 K (for B — 0). Therefore the experi- 
mentally observed increase of the escape rate with B is 
indeed due to the discussed mechanism of i?-enhanced 
tunneling. The smaller experimental values of the rel- 
ative escape rate W{B)/WiO) for T = 0.26 K can be 
understood by noticing that the activation rate is close 
to the tunneling rate for such T, and since it presum- 
ably only weakly depends on B , the overall slope of 
ln[W{B)/W{0)] should be smaller than that of the theo- 
retical curve which ignores activation (approximately, by 
a factor of 2, if we ignore the dependence of the activation 
rate on B). 



A. The prefactor 



The dependence of the potential U (z) ( |3q ) on n gives 
rise to the density dependence of the tunneling rate 
W{B) even for B = Q. We calculated the exponent 
and the prefactor in W{0) by matching the WKB wave 
function under the barrier for I/7 <C z <C L with 
the tail of the non-WKB intrawell solution (here, L = 
h^"/"^ /2m\e£±\ is the characteristic barrier width). In the 
spirit of the logarithmic perturbation theory (LPT) [ p6[ , 
the wave function of the ground state inside the well and 
not too far from it can be sought in the form 



4'g{z) — const X zexp[— A(z) 



(40) 



[we explicitly take into account that the function ipg{z) 
has a zero in the ground state]. 

The function dA/dz satisfies a Riccati equation. It can 
be solved near the well (z ^ L) by considering the last 
two terms in the potential U{z) (^) as a perturbation. 
For small z, the major correction comes from the term 
(X £±. To the first order in £j_, 



^(z)«7z(l 



(41) 



In obtaining this expression we took into account the cor- 
rection to the ground state energy 6Eg — — 3|e£^|/27. 
This correction can be obtained from the condition that 
the linear in £± term in dA/dz remain finite for z ^ 0. 

The correction to A (^Tj) is small for z small com- 
pared to the barrier width L. We note that the exponent 
A{z) has an overall functional form which differs from 
that of the commonly used [^ variational wave function 
tpiz) oc zexp(— 7z), with 7 being the variational param- 
eter. 
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FIG. 9. The rate of electron tunneling from the helium sur- 
face M^(0) for B = as a function of the electron density. The 
dots show the experimental data [^] . The pulling field £±_ was 
calculated from the geometry of the experimental cell, the 
applied voltage, and the density. 

The expression for A (^Tj) matches the small-z/L ex- 
pansion of the action S of the WKB wave function under 
the barrier for L z 7^^. This allowed us to find the 
prefactor in the WKB wave function and in the tunnehng 
rate. The resulting tunneling rate is shown in Fig. ^. It 
fully agrees with the experiment (see also Ref. pOt). 

VIII. CONCLUSIONS 

It follows from the results of the present paper that 
interelectron momentum exchange in correlated 2D elec- 
tron systems leads to an exponentially strong change of 
the rate of tunneling decay in a magnetic field parallel to 
the electron layer. The mechanism is dynamical by na- 
ture. It depends on the interrelation between the char- 
acteristic momentum exchange rate and the reciprocal 
duration of tunneling in imaginary time 1/t/. 

For low temperatures, where escape occurs via tunnel- 
ing from the ground state, the tunneling rate is affected 
primarily by high-frequency in-plane electron vibrations. 
In their turn, these vibrations are determined by short- 
range order in the 2DES. Their frequencies are of the 
order of the characteristic zone-boundary frequency of 
the Wigner crystal Wp. If ^ the effect of the 

magnetic field on tunneling is nearly completely compen- 
sated in the case where the width of the tunneling barrier 
is small compared to the interelectron distance. 

At higher temperatures, the magnetic field may in fact 
enhance rather than suppress the rate of tunneling de- 
cay. The overall escape rate as a function of B and T is 
expected to display a number of other unusual features. 
These include switching from activated escape to tun- 
neling and vice versa, and switching between tunneling 
from the ground and excited states. These switchings 
have been analyzed for simple but realistic models of the 
tunneling barrier. 



Our results on the field dependence of the tunneling 
rate and its evolution with temperature are in full qual- 
itative and quantitative agreement with the existing ex- 
perimental data on tunneling from a strongly correlated 
2DES on helium with no adjustable parameters. 

The results also apply to 2DES in semiconductor het- 
erostructures. For correlated systems in semiconductors, 
tunneling has been investigated mostly for the magnetic 
field B perpendicular or nearly perpendicular to the elec- 
tron layer, cf. The data on tunneling in a field par- 
allel to the layer refer to high density 2DESs |^ , where 
correlation effects are small. 

The effect of a parallel magnetic field is most pro- 
nounced in systems with shallow and broad barriers U{z). 
For example, in a GaAlAs structure with a square barrier 
of width L ~ 0.1 /xm and height 7^/2to = 0.02 eV, for 
the electron density n = 1.5 x 10^° cm^^ and i? = 1.2 T 
we have ujpTQ « 0.6 and uJcTq ~ 1 (tq = mL/j is the tun- 
neling duration for n = B = 0). The results of Sec. VI C 
for square barriers (with account taken of the correlation- 
hole correction) show that the interelectron momentum 
exchange should significantly modify the tunneling rate 
in this parameter range, provided the 2DES is correlated 
p9[ . We therefore expect that tunneling experiments on 
low-density 2DESs in parallel fields will reveal electron 
correlations not imposed by the magnetic field, give in- 
sight into electron dynamics, and possibly even reveal 
a transition from an electron fiuid to a pinned Wigner 
crystal with decreasing n. 
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